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ABSTRACT 

Early-type galaxy velocity dispersions and luminosities are correlated. The correlation estimated in 
local samples (< 100 Mpc) differs from that measured more recently in the SDSS. This is true even 
when systematics in the SDSS photometric and spectroscopic parameters have been accounted- for. 
We show that this is also true for the ENEAR sample if galaxy luminosities are estimated using 
distances which have been corrected for peculiar motions. We then show that, because the estimate 
of the 'true' distance is derived from a correlation with velocity dispersion, in this case the £)„ — a 
relation, using it in the cr — L relation leads to an artificially tight relation with a biased slope. Making 
no correction for peculiar velocities results in a cr — L relation which is very similar to that of the 
SDSS, although with larger scatter. We also measure the a — L correlation in a mock ENEAR catalog, 
in which the underlying galaxy sample has the same a — L correlation as seen in the SDSS. The mock 
catalog produces the same I?„ — a relation as the data, the same biased slope when Z?„ — a distances 
are used to estimate luminosities, and good agreement with the input a — L relation when redshift is 
used as the distance indicator. This provides further evidence that the true a — L relation of ENEAR 
galaxies is indeed very similar to that of SDSS early-types. Our results suggest that local a — L 
relations which are based on Fundamental Plane distances should also be re-evaluated. Our findings 
also have important implications for black hole demographics; the best direct estimates of the masses 
of supermassive black holes come from local galaxies, so estimates of the black hole mass function are 
more safely made by working with the M, — a correlation than with M, — L. 

Subject headings: galaxies: elliptical — galaxies: fundamental parameters — galaxies: photometry — 
galaxies: spectroscopy — black hole physics 



1. INTRODUCTION 

The luminosities and velocity dispersions of early-type 
galaxies are strongly correlated: their logarithms follow 
an approximately linear relation (e.g. Faber & Jack- 
son 1976). There is a long and complicated history of 
what the slope of this relation is, mainly due to the 
fact that, if there is intrinsic scatter around this cor- 
relation, then there are at least three interesting slopes: 
that obtained from fitting the mean luminosity at each 
velocity dispersion, (logLj logcr) = a^io- log a -|- 6L|cr, the 
mean velocity dispersion as a function of luminosity, 
(logCTjlogL) — a„\L\ogL 4- and the slope of the 

principal axis of the joint distribution of L and a (e.g. 
Lynden-Bell et al. 1988; SagUa et al. 2001). The first of 
these allows one to use a to predict L, whereas the sec- 
ond must be used if one wishes to predict a from L. Until 
recently, authors were not careful to distinguish between 
these cases. 

Because L depends on the distance to the source, 
whereas a does not, it is straightforward to estimate 
(logallogL) from flux hmited samples; naive estimates 
of the other two correlations are compromised by se- 
lection effects. Most studies of this correlation, based 
on local samples, agree that the typical velocity disper- 
sion at fixed luminosity, (logcrl logL), scales as L^/**: 
0'a\L ~ 1/4- The common unfortunate abuse of jargon 
is to say that luminosity scales as the fourth power of 
velocity dispersion; in fact, a]^\„ is considerably shal- 
lower. In what follows, we will be almost exclusively 
concerned with how well L can be used as a predictor for 
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cr: (logCT|logL). 

Prior to the year 2000, studies of the a — L correlation 
were restricted to local samples (from within about 50 
Mpc) containing ~ 100 objects. The SDSS provided a 
sample of early-type galaxies that was larger by about 
two orders of magnitude (Bernardi et al. 2003a): ^ 9000 
objects drawn from a volume which extended to consider- 
ably larger distances (median redshift z 0.1). However, 
the a — L relation in this sample. 



logioCr|Mr 



2.203- 0.102 (A/^ + 21), (1) 



SDSS-B03 

with an intrinsic scatter of about 0.07 dex is inconsis- 
tent with that found in local samples. For example, at 
log]^Q(cr/km s~^) > 2.4, the relation which Forbes & Pon- 
man (1999) found best fits the Pruniel & Simien (1996) 
sample of 236 local early-type galaxies differs from the 
SDSS fit by more than 0.05 dex. Expressed in terms 
of absolute magnitudes, the fits differ by more than 
0.5 mags at < —23. This is substantially larger than 
expected given the measurement errors. What causes 
this difference? 

Because the SDSS luminosities and velocity dispersions 
are obtained from an automated pipeline, (i.e., light pro- 
files and spectra were not inspected individually) , it may 
be that, say, the SDSS measurements differ systemati- 
cally from those estimated for local objects in the liter- 
ature. Indeed, it is known that the SDSS photometric 
reductions underestimate the luminosities of objects in 
crowded fields (Mandelbaum et al. 2005; Bernardi et 
al. 2006a; Hyde et al. 2006; Lauer et al. 2006). See 
the Appendix for a more detailed discussion. The Ap- 
pendix also shows that the SDSS also slightly overesti- 
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mates the velocity dispersions at small a. In what fol- 
lows, we will refer to the sample in which these system- 
atics have been accounted-for as the SDSS-B06 sample 
(see the Appendix). In the SDSS-B06 sample, 



log(T|A/r 



SDSS-B06 



2.190 - 0.100 + 21) (2) 



with intrinsic scatter of 0.07 dex; it happens that this is 
not very different from the relation obtained by Bernardi 
et al. (2003b). 

A common problem from which all local samples suf- 
fer is that, while the apparent magnitude of an object 
can be measured quite accurately, the absolute magni- 
tude is more difficult because it depends on the distance 
to the galaxy. The true distance is difficult to measure 
because the redshift, which is usually well-determined, 
is a combination of the distance to the galaxy and the 
component of its peculiar velocity which is directed along 
the line of sight to the observer: cz = Hd + vios- Typi- 
cal velocities are expected to be of order a few hundred 



km 



so the redshift is a reliable distance indicator 



only beyond about 100/i~^Mpc. Most of the objects in 
the SDSS lie well beyond this distance, whereas all lo- 
cal samples are shallower. Thus, a legitimate concern is 
whether uncertainties in estimating the true distance are 
driving the difference between the SDSS measurement 
and those which are based on more local samples. 

To address such concerns, we have studied the a — L 
relation in the definitive sample of nearby early-type 
galaxies — that assembled in the ENEAR database (da 
Costa et al. 2000; Bernardi et al. 2002a; Alonso et al. 
2003; Wegner et al. 2003). ENEAR contains about 1000 
objects out to 7000 km s"-^, for which measured redshifts 
and estimated distances are available. Section 2 shows 
how the (log (T I log L) relation changes if one uses the 
redshift rather than the estimated distance when com- 
puting L. It highlights the fact that, if a played a role 
in determining the distance, e.g. if the distance comes 
from a Fundamental Plane or D„ — a analysis, then the 
a — L relation may be biased. In Section 3 we discuss 
our results. 

2. THE LOCAL a - L RELATION 

The first part of this section compares various deter- 
minations of the local a — L relation. The second part 
compares these determinations with that from the SDSS- 
B06. 

2.1. ENEAR 

The ENEAR sample (da Costa et al. 2000) is approx- 
imately magnitude limited to 14.5 in the B band. The 
catalog provides new apparent magnitudes, dn measure- 
ments, redshifts, and velocity dispersions for about 1000 
early-type galaxies distributed over the whole sky out to 
about 7000 km s"^ 

In principle, estimating {(j\L) in the ENEAR sample is 
straightforward because the ENEAR team has also pub- 
lished a measurement of the Z?„ — a correlation: 
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Fig. 1. — Joint distribution of L and a in the ENEAR sam- 
ple, when L has been estimated using a distance which has been 
corrected for peculiar motions. Short dashed line shows {a\L^). 
Using the redshift as distance indicator instead leads to the {a\Lz) 
relation shown as the solid line. The symbols and lines show the 
result of using total rather than bulge luminosities for ENEAR. 

(Bernardi et al. 2002b), with a scatter of 0.2/ln(10) dex. 
Here, dn is that angular scale within which the average 
surface brightness in the Rc band is 19.25 mag/arcsec^. 
For reasons we describe below, we compute luminosities 
for the ENEAR galaxies based on two different distance 
estimates. One uses this D„— cr relation to estimate 'true' 
distances, and the other uses the redshift as a distance 
indicator (i.e., this second distance estimate ignores pe- 
culiar velocities). We will refer to the associated lumi- 
nosities as Ld and Lz- In addition, for ease of comparison 
with the SDSS, we shift aU magnitudes from ENEAR Rc 
to SDSS r = 0.24 + Rc. 

The joint distribution of a and Ld which results from 
using the Dn—cr distance estimate (filled circles) is shown 
in Figure [H Here we only show galaxies which have both 
spectroscopy and photometry observed by the ENEAR 
team (i.e. we did not include measurements listed in 
the ENEAR catalog which were compiled from previous 
work). In addition, we selected galaxies with disk-to- 
bulge ratio smaller than 0.5 and used the total magnitude 
fisted by Alonso et al. (2003) (using the bulge luminosity 
the results only slightly change). Dashed line shows 



log 



10 



/ ENEARd 



2.159 - 0.130 (M^ + 21) (4) 



which best-fits the sample. Comparison with equa- 
tion ^ shows that this relation is considerably steeper 
than in the SDSS. 

The observed scatter around this relation is ^ 0.08 dex. 
This is about 0.01 dex larger than what SDSS reports 
as intrinsic scatter. This is remarkable, because ENEAR 
ought to be carrying distance errors of about twenty per- 
cent. These errors ought to translate into increased scat- 
ter of 0.43 mags along the x-axis of Figure [T] If not 
accounted-for (and the fit above does not), this should 
have artificially decreased the slope and increased the 
scatter of the {a\L) correlation. E.g., adding this scatter 
to the SDSS numbers would have decreased the magni- 
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Fig. 2. — Joint distribution of L and c in a mock catalog in which the underlying galaxy population has been constructed to mimic the 
SDSS-B06 sample, and to which the ENEAR, magnitude and redshift selection cuts have been applied. The two panels differ in how the 
luminosities were estimated: the panel on the left uses the observed redshift to define a distance, whereas the panel on the right makes 
a correction for peculiar velocities which is based on the D„ — cr relation. In both panels, solid line shows the input {fT\L) relation, and 
dashed line shows the relation defined by the points in the panel. 



tude of the slope by about fifteen percent, and the rms 
scatter around the new relation would have been about 
fifteen percent larger. 

2.2. To remove Wpoc or not to remove tipec ? 

Since the I?„ — a estimate of the distance is only 
good to 20 percent, one might well wonder if sam- 
ples like ENEAR reach distances at which the red- 
shift itself provides a more reliable estimate of the 
true distance than relations like I?„ — cr or the Fun- 
damental Plane. This will happen on scales r where 
(upcc/300 km s"^)/0.2 < (r/S/i-^Mpc). This suggests 
that, beyond about 50/i^^Mpc, the redshift may actu- 
ally be a better estimate of the true distance even if 
typical line-of-sight peculiar velocities were as high as 
1000 km s"^ 

Since ENEAR straddles this distance regime, we have 
studied what happens to the a — L relation when the 
redshift itself is used as an estimate of the distance — 
i.e., no correction is made for the peculiar velocity. The 
inferred relation, when the total luminosity is used. 



a\M) 



= 2.184- 0.104 + 21), (5) 

/ ENEARz 

is shown as the dashed line in Figure [T] The relation 
obtained when the total luminosity is scaled to the bulge 
luminosity 

flogin cr|M\ = 2.198 - 0.102 (AU + 21), 

\ / ENEARz-bulgc 

is shown as the dotted line. The two relations are similar. 
Equation [5] is considerably shallower than {a\Ld) (equa- 
tionHI, with substantially larger scatter (0.12 dcx). 

Comparison of the relation based on redshift, a — 
say, with the one based on distances, a — Ld, shows that 



the scatter appears to decrease considerably when us- 
ing the distances rather than the redshifts. While it is 
tempting to conclude that this is signalling that the dis- 
tance indicator is accurate, this is not the full story. Af- 
ter all, the Dn — a distance estimate is actually a func- 
tion of a: distance cx a^-"^ jdn. As a result, a appears 
in both axes of Figure [U with the x-axis proportional 
to m — 5 log]^o('''^'^/'^n)- Therefore, scatter in the re- 
lation will be correlated along a line which has slope 
y = — 5(1.2x). A slope of —1/6 is not far-off the one actu- 
ally observed. Hence, it may be that the relation shown 
in Figure [Hand quantified by equation (j4]) is both biased 
and artificially tight. Note that similar concerns would 
also apply to the use of Fundamental Plane-based dis- 
tance estimates (such as those used by Pruniel & Simien 
1996), since the coefficient of the velocity dispersion in 
the direct fit is also ^ 1.2 (e.g. J0rgensen et al. 1996; 
Bernardi et al. 2003c). This discussion provides another 
reason why the simple procedure of using the redshift 
as the distance indicator might actually be the preferred 
one. 

2.3. Results from a mock ENEAR-SDSS catalog 

To better understand the effect of using a cr-based dis- 
tance estimate on the a — L relation, we constructed a 
mock catalog of the ENEAR sample. This was done by 
assuming that the joint distribution in luminosity, size 
and velocity dispersion for ENEAR galaxies is the same 
as for SDSS-B06 galaxies, so we could simply follow the 
steps described by Bernardi et al. (2003b). We then 
added Gaussian noise with dispersion equal to the EN- 
EAR observational estimates, and applied the ENEAR 
apparent magnitude and redshift cuts. 

Since this method allows us to generate both L and a 
half-light radius Re, the assumption that the light pro- 
file is deVaucoleur allows us to compute a value of d„ 
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Fig. 3. — Joint distribution of L and a in the ENEAR sample when the luminosity is computed using the redshift distance; i.e., the 
distance has not been corrected for peculiar motions. Short dashed lines show the {a\Lz) relation defined by these points, and solid line 
shows the relation defined by the SDSS-B06 sample. Dotted line shows the result of using bulge rather than total luminosities for ENEAR. 



for each mock galaxy. We then fit for the D„ — a rela- 
tfon, finding that it had slope 1.15, zero-point 1.4 and 
rms scatter 0.1 dex. These values are rather similar to 
equation ([3]), suggesting that our mock catalog is actu- 
ally rather realistic. 

Figure [5] compares the a — Ld and a — Lz relations in 
the mock catalog; the relation based on the distance indi- 
cator is clearly steeper and tighter than the one based on 
the redshift, just as in the ENEAR sample itself. More- 
over, notice that the slope of {<j\Ld) relation is clearly 
steeper than the input slope, whereas {a\Lz) is in rather 
good agreement with the input relation. This suggests 
that equation ([5]) is closer to the true relation for EN- 
EAR galaxies than is equation Q. This bias does not 
depend strongly on the range of apparent magnitudes 
in the sample. It does depend on the intrinsic distri- 
bution of the absolute magnitudes. When the absolute 
magnitude range is large then the slope of the a ~ Ld re- 
lation is closely related to that of the D„ — a relation (see 
discussion in previous subsection) and has little scatter. 
However, when the luminosity function is narrow (i.e. 
sharply peaked) then the a — Ld relation shows larger 
scatter: while the slope is still biased it is less easily re- 
lated to that of the £>„ — a relation. 

2.4. Comparison with SDSS-B06 

The analysis above suggests that the a — Lz relation 
in ENEAR, and, by extension, other local samples, are 
likely to be more reliable than the local a — Ld rela- 
tions. Figure [3] shows the joint distribution of Lz and 
a in the ENEAR sample when the luminosity is com- 
puted using the redshift distance; i.e., the distance has 
not been corrected for peculiar motions. The dashed line 



shows equation ([5]), and the dotted line shows the result 
of using bulge rather than total luminosities for ENEAR. 

These fits should be compared to the solid line, which 
shows the SDSS-B06 relation (equation [2]), for which the 
redshift is an excellent indicator of the true distance. 
Note that the ENEAR and SDSS-B06 samples are in 
excellent agreement, suggesting that correcting for pecu- 
liar motions in local samples can lead to serious biases in 
correlations which involve cr. 

3. DISCUSSION 

We showed that the ENEAR a — Ld relation, where 
Ld indicates that the luminosity is based on a distance 
estimate which has been corrected for peculiar motions, 
is almost certainly biased. This is because the distance 
indicator used to estimate Ld depends on a. We also 
showed that if the redshift is used as the distance indica- 
tor (i.e., no correction for peculiar motions is made) then 
the resulting a — Lz relation, while noisier, is less likely 
to be biased. The ENEAR a — Lz relation is in excellent 
agreement with that measured in the SDSS, where ne- 
glecting peculiar velocities is an excellent approximation. 
Thus, it appears that the discrepancy between the SDSS 
a — L relation and that in local samples is due to the use 
of a distance indicator which correlates with a. 

Our results are of particular interest for the problem 
of estimating black hole abundances from the observed 
distribution of luminosities or velocity dispersions (e.g. 
Yu & Tremainc 2002; Shankar et al. 2004; Tundo et al. 
2006). The first method requires {M,\L) whereas the sec- 
ond requires (Af, jcr). The two approaches will only give 
the same estimate of 4){M,) if the {a\L) relation of the 
black hole sample is the same as that of the sample from 
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which the luminosity function (t>{L) and velocity disper- 
sion function 4>{a) are drawn. Current black hole samples 
are relatively local, so their a — L correlations are very 
different from that of the SDSS. Bernardi et al. (2006b) 
suggests that current black hole samples are biased to- 
wards objects with abnormally large velocity dispersions 
for their luminosities. If this is a selection rather than 
physical effect, then the M, — a and M, — L relations 
currently in the literature are also biased from their true 
values. Bernardi et al. find that the bias in the (M, |cr) 
relation is likely to be small, whereas the (M, \L) relation 
is biased towards predicting more massive black holes for 
a given luminosity. Therefore, the estimate based on the 
velocity function (e.g. that of Sheth et al. 2003) is to be 
preferred. 

The (cr|L) relations show evidence for a flattening at 
large L, a fact we do not use here, but which is relevant 
to studies of BCGs (e.g. Bernardi et al. 2006a) and may 
be relevant to studies of black hole demographics. 

I would like to thank Ravi Sheth for discussions on 
how to generate the ENEAR mock catalog and Moniquc 
AUer for comparisons of velocity dispersion measure- 
ments. This work is partially supported by NASA grant 
LTSA-NNG06GC19G, and by grants 10199 and 10488 
from the Space Telescope Science Institute, which is op- 



erated by AURA, Inc., under NASA contract NAS 5- 
26555. 

Funding for the SDSS and SDSS-II has been pro- 
vided by the Alfred P. Sloan Foundation, the Partici- 
pating Institutions, the NSF, the US DOE, NASA, the 
Japanese Monbukagakusho, the Max Planck Society and 
the Higher Education Funding Council for England. The 
SDSS website is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions: 
The American Museum of Natural History, Astrophys- 
ical Institute Postdam, the University of Basel, Cam- 
bridge University, Case Western Reserve University, the 
University of Chicago, Drexel University, Fermilab, the 
Institute for Advanced Study, the Japan Participation 
Group, the Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Parti- 
cle Astrophysics and Cosmology, the Korean Scientist 
Group, the Chinese Academy of Sciences (LAMOST), 
Los Alamos National Laboratory, the Max Planck Insti- 
tute for Astronomy (MPI-A), the Max Planck Institute 
for Astrophysics (MPA), New Mexico State University, 
the Ohio State University, the University of Pittsburgh, 
the University of Portsmouth, Princeton University, the 
U.S. Naval Observatory, and the University of Washing- 
ton. 



REFERENCES 



Adelman-McCarthy, J. et al. 2006, ApJS, 162, 38 

Alonso, M. v., Bernardi, M., da Costa, L. N., Wegner, C, Willmer, 

C. N. A., Pellegrini, P. S., & Maia, M. A. G. 2003, A J, 125, 2307 
Bernardi, M., Alonso, M. V., da Costa, L. N., Willmer, C. N. A., 

Wegner, C, Pellegrini, P. S., Rite, C, & Maia, M. A. G. 2002a, 

AJ, 123, 2990 

Bernardi, M., Alonso, M. V., da Costa, L. N., Willmer, C. N. A., 
Wegner, G., Pellegrini, P. S., Rite, C, & Maia, M. A. G. 2002b, 
AJ, 123, 2159 

Bernardi M., Sheth R. K., Annis J., et al. 2003a, AJ, 125, 1817 
Bernardi, M., Sheth, R. K., Annis J. et al. 2003b, AJ, 125, 1849 
Bernardi, M., Sheth, R. K., Annis J. et al. 2003c, AJ, 125. 1866 
Bernardi, M., Hyde, J. B.. Sheth. R. K.. Miller, C. J., Nichol. R. 

C. 2006a, AJ, in press (astro-ph/0607117) 
Bernardi, M., Sheth. R. K.. lundo, t^., h Hyde, J. B. 2006b, ApJ, 

submitted (astro-ph/0609300 ) 
da Costa, L. N., hSernardi, M., Alonso, M. V., Wegner, G., Willmer, 

C. N. A., Pellegrini, P. S., Rite, C, & Maia, M. A. G. 2000, AJ, 

120, 95 

Desroches, L.. Ouataert. E.. Ma, C, & West, A. A. 2006, MNRAS, 

submitted (astro-ph/0608474) 
Faber, S. M. & Jackson, K. 1976, ApJ, 204, 668 
Forbes, D. A., & Ponman, T. J. 1999, MNRAS, 309, 623 



Hyde, J. B., Bernardi, M., Sheth, R. K. et al. 2006, AJ, submitted 
j0rgensen, I., Franx, M., & Kjasrgaard, P. 1996. MNRAS. 280, 167 
Lauer, T. R., et al. 2006, ApJ, submitted ( astro-ph/06_06739 ) 
Lynden-Bell, D., Faber, S. M., Burstein, D., uavies, it. Li., uressler, 

A., Terlevich, R. J., & Wegner, G. 1988, ApJ, 326, 19 
MitcheU, J. L., Keeton, C. R., Frieman, J. A., & Sheth, R. K. 2005, 

ApJ, 622, 81 
Prugniel, P., & Simien, F. 1996, A&A, 309, 749 
Sagha, R. P., CoUess, M., Burstein, D., Davies, R. L., McMahan, 

R. K., & Wegner, G. 2001, MNRAS, 324, 389 
Shankar, F., Salucci, P., Granato, G. L., De Zotti, G., & Danese 

L. 2004, MNRAS, 354, 1020 
Sheth, R. K., Bernardi, M., Schechter, P. L., et al., 2003, ApJ, 594, 

225 

Tundo, E., Bernardi, M., Hvde. J. B.. Sheth, R. K., & Pizzella, A. 

2006, ApJ, submitted ( astro-ph/0609297 1 
Wegner, G., Bernardi, M., Willmer, C. N. A., da Costa, L. N., 

Alonso, M. v., Pellegrini, P. S., & Maia, M. A. G. 2003, AJ, 126, 

2268 

Yu, Q., & Tremainc, S. 2002, MNRAS, 335, 965 



APPENDIX 

COMPARISON OF THE a - L RELATION IN VARIOUS SDSS DATA RELEASES 

Bernardi et al. (2003b) report the first SDSS-based a — L relation, from a sample of about 9000 objects classified as 
being early- types on the basis of imaging (concentrated light profiles) and spectroscopy (weak or no emission lines) : 
hereafter, we will refer to this sample as SDSS-B03. 



Changes to SDSS photometry 

Shortly after this estimate was published, the SDSS reported a problem with the point-spread-function in the SDSS 
photometry which led, on average, to an overestimate of the apparent brightness of extended objects and larger effective 
radii. This problem was fixed in subsequent SDSS data releases. However, more recent work has shown that there 
remain a problem with the SDSS photometric reductions which is most severe for bright objects in crowded fields 
(the SDSS pipeline overestimates the sky level; see Mandelbaum et al. 2005; Bernardi et al. 2006a; Lauer et al. 
2006). Re-analysis of such objects (Bernardi et al. 2006a; Hyde et al. 2006) suggests that the SDSS-DR5 photometric 
reductions underestimate the true apparent brightness by 0.1 mags on average, and by up to 0.5 mags for bright objects 
in crowded fields (by chance the combination of the two problems make the magnitudes in Bernardi et al. 2003a to 
be similar to the recomputed values by Hyde et al.). 
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These changes in the photometric reductions affect the estimated L in the a — L relation. They also have a small 
effect on <j, because the velocity dispersions are scaled (to account for aperture effects) to a fraction (typically 1/8) 
of the half light radius, and the half light radii are about 15% smaller than in Bernardi et al. (2003a). This aperture 
correction is small, so the change to the a — L relation is driven by the change to L. 

Changes to SDSS spectroscopy 

In addition, the SDSS spectroscopic reductions app ear to have changed between Bernardi et al. (2003a) and DR5. 
The left-hand panel of The top left panel in Figure [Al] shows that cr in the SDSS-DR5 for the SDSS-B03 sample do not 
match the values used by Bernardi et al. (2003a). The difference is small but systematic, with SDSS-DR5 being larger 
than SDSS-B03 at small cr. Since there is considerably more overlap between measurements from the literature (from 
the HyperLeda database) with SDSS-DR5 than with SDSS-B03, we have checked if the SDSS-DR5 velocity dispersion 
values are biased at low a. The top right hand panel of Figure lAll shows that SDSS-DR5 is biased towards larger a 
at (7 < 150 km s~^; The bias is similar to that seen on the top left panel. 

The SDSS-DR5 values are based on averaging the measurements from the Direct- Fitting and Fourier-fitting methods 




Fig. Al. — Comparison of cr values from Bernardi et al. (2003a) (left panels) and HyperLeda (right panels) with those from SDSS-DR5 
(top) and with the values used in this paper (bottom). Top left: At fixed SDSS-B03, the SDSS-DR5 values are systematically higher when 
<j < 150 km s— 1. Bottom left: The same, but after correcting for this bias by re-analyzing the spectra. Top right: At small cr, SDSS-DR5 
reports larger values than HyperLeda; the offset is similar to that seen in the left-hand panel; in this case, we know that SDSS-DR5 is 
biased high. HyperLeda reports substantially larger velocity dispersions at ct > 250 km . The following figures suggest that, in this 
case, SDSS is more reliable than HyperLeda. Bottom right: The same, but after correcting for this bias by re-analyzing the spectra. 
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Fig. A2. — Comparison of the SDSS observed spectrum of an object in common with the hterature (jagged hne) with templates broadened 
by the velocity dispersions reported by SDSS-DR5 (288 km s~^) and HyperLeda (436 km s~^) (red and green lines respectively). Top and 
bottom panels show two sections of the spectrum. The smaller velocity dispersion is clearly a better fit. 




Fig. A3. — Cross correlation peak for the object shown in the previous Figure. The SDSS-DR5 value is clearly a better description of this 
broadening function (red lines was computed using the SDSS velocity dispersion while green is from the value in the literature). Panel also 
shows estimates of the velocity dispersion obtained from the Direct Fitting, Fourier Fitting and cross-correlation methods: the SDSS-DR5 
pipeline returns the average of the first two methods. 



(Bernardi et al. 2003a). We have run simulations similar to those in Bernardi at al. (2003a) and found that the 
discrepancy results from the fact that the Fourier-fitting method is now biased 15% level at low sigma (~ 100kms~^), 
whereas the other method is not. The bottom panels of Figures lAll shows that using only the Direct- Fitting method 
improves the agreement between SDSS-DR5 and SDSS-B03 as well as between SDSS-DR5 and HyperLeda. 

The right panels in Figure lAll show that HyperLeda reports substantially larger velocity dispersions at cr > 
250 km s~^. Figures [A2I and |A3I suggest that, in this case, SDSS is more reliable than HyperLeda (these figures 
are a representative example; we have checked all galaxies with velocity dispersion larger than 250 km s~^). This is 
important, because the objects with the largest a are expected to host the most massive black holes. If the SDSS 
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Fig. A4. — Comparison o{ a — L relation in various SDSS reductions. Panel on the left shows the relations themselves, and panel on 
the right shows the relations expressed as residuals from the fit given in Equation ((2J. Although this is not the main point of the current 
study, note that the relation clearly flattens at large magnitudes. 

determinations were indeed biased low, then determinations of black hole abundances which are based on the SDSS 
velocity function (Sheth et al. 2003) would underestimate the abundances of supcrmassive black holes (e.g. McClure 
& Dunlop 2004; Shankar et al. 2004; Tundo et al. 2006). Our analysis shows that this is not a concern. 

The SDSS-B06 sample and the a — L relation 

Because of the problems with the SDSS photometry and spectroscopy, we have chosen to use our new estimate of the 
a — L relation. We selected the SDSS-B03 sample of early-type galaxies. For the photometry, we use the reductions 
from Hyde et al. (2006). For the velocity dispersions, we use the Direct- Fitting method described above, but do not 
average it with the Fourier-fitting method (as done for the values in the SDSS database). In the main text, we refer 
to the catalog which results as the SDSS-B06 sample^ 

The main text studies the a — L relation. Figure [X4l illustrates the effect of these various SDSS reductions on this 
relation. In all cases, straight lines show single power-law fits to (cli), and jagged lines show the median cr as a function 
of L for a few bins in L. To highlight the differences, the y-axis in the left hand panel shows a — 2.19 — 0.100(Afr + 21), 
rather than a itself, as a function of L. The triple-dot dashed lines show the relation reported by Bernardi et al. (2003b), 
and the dotted fines show fits based on quantities output by SDSS-DR5 (corrected for aperture effects to Re/8)- The 
short dashed lines are based on SDSS-DR5 photometry, but use the SDSS-B06 dispersions (to approximately correct 
for the small systematic bias in SDSS-DR5 us), and solid lines show the result of recomputing the photometry as well 
as the velocity dispersions (the SDSS-B06 sample). 

The relations from SDSS-B03 and SDSS-B06 are rather similar (the various adjustments to the SDSS-B03 photometry 
approximately cancel). Notice that all these relations show evidence for a flattening at large L, a fact we do not use 
in the main text, but which is relevant to studies of BCGs (e.g. Bernardi et al. 2006a) and may be relevant to studies 
of black hole demographics. 

Object by object comparison with ENEAR 

The main text compares the a — L relation in the SDSS-B06 and ENEAR samples. In the remainder of this section 
we provide an object by object comparison of the two sets of photometric and spectroscopic reductions. We began 
by matching the two catalogs: only about 50 objects with ENEAR imaging and/or spectroscopy were also in the 
SDSS-DR5 sample. This surprisingly small number is because the ENEAR objects are, in general, very bright and 
very large compared to the vast majority of objects targeted for SDSS spectroscopy. Of these, 30 have ENEAR velocity 
dispersions and 29 have ENEAR photometry, but only 15 have both. 

Figure [XSl compares the ENEAR and the Direct- Fitting estimates of the velocity dispersions of the 30 objects which 
have ENEAR spectroscopy available. The estimates are generally in good agreement, suggesting that SDSS-B06 
velocity dispersions are unlikely to be systematically biased with respect to local samples. Therefore, the discrepancy 
in the a — L relation is almost certainly due to systematic errors in the luminosity. 

Such errors could arise from systematic differences in the photometry itself, or in the conversion from apparent 
to absolute magnitude. To test this. Figure IA5I compares the ENEAR and Hyde et al. estimates of the apparent 
magnitudes for 29 galaxies which have ENEAR imaging available. Whereas Hyde et al. reports a magnitude associated 
with fitting a single deVaucoleur profile to the image, ENEAR reports both total (filled circles) and bulge (open 
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Fig. A5. — Left: Comparison of ENEAR and Direct- Fitting estimates of the velocity dispersions of the 30 objects in common. There 
is no offset, and the magnitude of the scatter is consistent with being entirely due to observational errors. Right: Comparison of ENEAR 
and Hyde et al. estimates of the apparent magnitudes of the 29 objects in common. Whereas Hyde et al. reports a magnitude associated 
with fitting a single deVaucoleur profile to the image, ENEAR reports both bulge (open diamonds) and total (filled circles) magnitudes. 
On average, the Hyde et al. magnitude is similar to ENEAR total magnitude, but is 0.12 mags brighter than ENEAR bulge magnitude. 
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Fig. A6. — Joint distribution of ENEAR (open diamonds) and the new recomputed SDSS (filled circles) luminosity and velocity dispersion 
estimates for the 15 objects common to both surveys with available imaging and spectroscopy. Dashed and solid lines show {<j\L) for the 
full ENEAR and SDSS-B06 samples, when the luminosities are estimated from the redshifts (i.e. no correction for peculiar velocities is 
made). 



diamonds) magnitudes. The figure indicates that Hyde et al. magnitudes tend to be similar to the ENEAR total 
magnitudes, and about 0.12 mags brighter than ENEAR bulge magnitudes. Here we compare galaxies dominated by 
the bulg e com ponent (i.e. with disk-to-bulgc ratio less than 0.5). 

Figure ED shows the a — Lz relation for objects with both ENEAR imaging and spectroscopy in common with SDSS. 
It is in reasonably good agreement with the SDSS-B06 relation. 



